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Morphology and testing of the retinal fundus images since it has the fastest execution 
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Network images are performed for the classification and grading of retinal fundus 


images into normal and abnormal using Alyuda Neuro-Intelligence software. 
A patient database is created using MySQL to store the required details of the 
patient and a graphical user interface is developed for an efficient usage of 
the system. The execution time of the system is found to be 7-9 seconds and 
is tested on 270 retinal fundus images. The precision and accuracy of the 
algorithm is 92.5% and 93.9%, respectively. 
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1, INTRODUCTION 

Diabetes essentially happens when the pancreas fails to secrete sufficient insulin for metabolism. It 
is a long term condition that causes very high glucose levels. Diabetic Retinopathy (DR) is a leading eye 
related disorder affected by diabetes. Since, it is asymptomatic in the initial stage; an efficient screening will 
prevent blindness. Various Image processing algorithms are used for the detection of Microaneurysms and 
Hemorrhages, which are the early clinical signs of DR. The detection of the early DR lesions and 
classification will help in combating blindness across the world. 

It is highly important for diabetic patients to have regular eye check-ups. Current retinal 
examination systems, which are utilized for locating and reviewing diabetic retinopathy, incorporate 
Ophthalmoscopy (which may be indirect and direct), Fluorescein angiography, and Fundus photography. 
These methods of assessment and detection of DR are currently manual, expensive and requires trained 
ophthalmologists. There is significant need of ophthalmologists across the world, especially in the rural 
areas. A computer aided diagnosis (CAD) system would tackle the problem by increasing the mass screening 
of the early DR patients before consulting the ophthalmologists. This in turn reduces the time consumption 
and increases the efficiency of the diagnosis of the disease. 
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The retinal features are macula, blood vessels, fovea and optic disc (OD). Any variations in these 
features may result in different retinal abnormalities. DR is of two types, namely Non proliferative DR 
(NPDR) and proliferative DR (PDR). Many features such as Microaneurysms (MAs), Hemorrhages (HMs), 
Exudates (hard and cottonwool spots) and Intraretinal Microvascular Abnormalities (IRMA) exists in NPDR. 
The very first sign related to early DR are microaneurysms (MAs) associated with the local expansions of the 
retinal capillaries and that result in hemorrhages (HA) when ruptured. In PDR, very small blood vessels 
develop from the retinal surface. The retinal fundus image with different features are as shown in Figure 1. 





Figure 1. Abnormal Retinal Fundus Image 


This paper proposes a CAD system for the early detection of DR to detect and characterise the red 
lesions based on a Hybrid technique clubbed with neural networks, a GUI that is user friendly and a patient 
database created using MySQL. The paper is organised in this fashion. A concise study of the research linked 
to the proposed system is presented in Section 2. A proficient early DR prediction system is presented in 
Section 3. The comprehensive results and discussions are given in Section 4. Conclusions are presented in 
Section 5 that sums up the explored work. 


2. RELATED WORK 

Most early screening programs use non-mydriatic digital fundus cameras to capture the retinal 
images. Retinal specialists then examine these images for the presence of lesions that are indicative of DR. 
However, the weak link is that the number of experts are considerably less and evaluating thousands of 
images is practically difficult. Thus, automatic detection via application of image processing techniques 
reduces the burden of the human experts.MA and HA are the earliest signs of DR. Blood vessels and OD 
needs to be eliminated for successful detection of early DR. 

Blood vessels needs to be eliminated from the retinal images since it has similar intensity as that of 
microaneurysms. Annie et.al [1] proposed an algorithm for retinal blood vessel segmentation, which 
consisted of pre-processing stages, enhancement using fast discrete curvelet transform and multi structural 
element morphology to detect retinal blood vessels. Curvelet transform was then applied on the retinal fundus 
images. The limitations include the tedious calculation of the vascular parameters. Soares et.al [2] suggested 
a method for the vessel extraction using Gabor wavelet and Bayesian classifier. The accuracy was found to 
be 96%. The method produced segmentation by categorizing every image pixel (as vessel or nonvessel), 
centered on pixel's features. Usman et.al [3] suggested a wavelet based method for the enhancement of the 
vessel along with the segmentation to obtain the region of interest. 2-D Gabor wavelet to augment the less 
visible vessels were used. Sensitivity was found to be around 94%. Similar works were carried out by [4], [5] 
and [6]. Ben et.al [7] proposed a method based on morphological scale space. Line structuring element was 
rotated around the seed point to obtain the curvature of the vessels. It was performed on few images, which 
was one of the main limitation. Zhu et.al [8] proposed a supervised approach for identifying the retinal 
vessels based on extreme learning machine for pixel classification. Average accuracy was around 96%. 
Hybrid classifiers could be used to improve the performance. Vasanthi et.al [9] carried out similar works. 
Akhavan et.al [10] suggested a method for the identification of vessels using fuzzy segmentation. The 
method showed consistent performance for normal as well as abnormal images. 

Detection and removal of OD is also essential for the early detection of DR. Various methods have 
been proposed for the elimination of OD. Ana et al. [11] proposed a technique for the extraction of the 
vascular tree using a technique based on graph-cut approach. This information was used to locate the OD. 
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The sensitivity was found to be less, which was a drawback. Akyol et. al. [12] proposed a technique for the 
identification of optic disk that included 

image-processing, key-point extraction, texture-analysis, visual-dictionary, and classifier techniques. 
Accuracy of the method was found to be around 94%. The technique worked well for images with noise and 
lesions. Murti et.al [13] suggested a method for the recognition of OD using least square fitting algorithm. 
Enhancement was done using histogram equalization and thresholding. Accuracy was found to be around 
97.5%. Less sensitivity and specificity were its major drawback. Rama Krishnan et al. [14] did similar works. 

Microaneurysms and Haemorrhages were detected after the elimination of the Blood vessels and 
Optic disk. Different methods were adopted for the detection of these early DR lesions. Preeyaporn et al. [15] 
proposed a method that used a mix of HSV technique, identification of area along with eccentricity 
technique. A color bar was framed with the colors of arbitrarily identified MAs. The position of the MA was 
recognized if the dot in the image exists within the range of the target color. The accuracy was 93%. It was 
found to have less sensitivity. Rukmini et al. [16] proposed another technique in light of fractal examination 
for the location of MAs. Fractal measurement was figured utilizing Box-Counting system. If the fractal 
dimension of a given image was lower than the threshold value, it was considered as normal or vice-versa. 
The sensitivity and the specificity were 89.5% and 82.1%, which is the main limitation. 

Akara Sopharak et al. [17] proposed a system in view of morphological operations for finding MA. 
After pre-processing, coarse segmentation was performed to distinguish the MA. Weka information mining 
programming and Naive Bayesian arrangement were utilized. The estimations of accuracy and precision were 
83.34 % and 99.99 %, respectively and sensitivity was found to be less. More improved results were 
achieved by [18]. Istvan et.al [19] proposed a system for the recognition of microaneurysm taking into 
account nearby rotational cross-segment profile analysis. Peak recognition was done on the acquired cross- 
sectional profiles and afterwards factual measures were computed to get the MA. Naive Bayes (NB) classifier 
was utilized. The advantages were high sensitivity with low false positive rates. Marwan and Eswaran [20] 
proposed a system in which MAs and HAs were recognized for the early analysis of DR utilizing CLAHE 
and h-minima transform. Classifiers were not used, which was one of the major limitation. Niladri et al. [21] 
proposed a technique for the identification of MA, which comprised of division of the retinal image into four 
quadrants, followed by preprocessing in all the four quadrants, and later was concatenated. Execution time 
was one of the major limitations. 

The survey of the literature works demonstrates that noteworthy works has been done in the field of 
identification of DR. Anyway, a framework with an automatic and early DR diagnosis system 1.e., MA and 
HA as early DR signs is still not productive. The main goal of this work is development of a framework for 
the early recognition of DR. It expects to distinguish MAs, HAs from the normal images precisely, group and 
grade them utilizing the neural network system in conjunction with storage of the details in the patient 
database and a Graphical user interface(GUI) for the proficient and simple access of the framework. It 
focuses fundamentally on the development of a hybrid framework that can effectively locate the early DR 
features with an enhanced sensitivity, specificity, precision and accuracy. 


3. MATERIALS & METHODS 

For effective detection of red lesions, the contrast between red lesions and the retinal background 
should be high and contrast between the retinal background and bright lesions should be low. This results in 
efficiently reducing the false positives during segmentation of candidate red lesions. Local entropy 
thresholding (LET) technique, takes the spatial distribution of gray levels into consideration and efficiently 
distinguishes enhanced dark lesions and the background as it can preserve the structural details of an image. 
Combining the best features of Local entropy thresholding and morphology (facilitates detection of smaller 
vessels), Hybrid method is proposed for the early detection of DR as shown in Figure 2. The CAD system for 
the early diagnosis of DR is divided into five phases: Acquisition of the Images, Pre-processing, Post- 
processing which includes removal of features like Optic Disk, Blood Vessels, Exudates, Feature Extraction 
of MA and HA, Classification and storage of patient data in the database. 

The retinal images are taken from many retinal fundus databases like DIARETDBO, DIARETDBI1, 
and hospital databases of Prabha Eye Clinic and Narayana Nethralaya, Bengaluru. Then, green component of 
the image is taken since green color plane shows the best vessel contrast. Red Lesions appear brighter in the 
green plane. Optic Disk has features like constant size, high intensity and circular shape. Initially, the 
maximum value of the green component of the image and the median is calculated. Then, a circular mask is 
constructed at the center co-ordinates using (1). 
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Figure 2. Block Diagram of the Hybrid method 


(a-p)* + (b-q)*=r" (1) 


Where, r is the radius and (p, q) are the center co-ordinates. The circular mask is then overlaid on the retinal 
image to get the OD masked image. Thus, Optic Disk is removed from the image. Background removal 
mainly aims at removing the background variations in the image. The foreground features are more 
prominent. Top and Bottom Hat Transform is applied to the OD masked image (foa). Top-hat operation is 
subtracting the result of performing a morphological opening on the OD masked image foa from the OD 
masked image (foa) using structuring element (SE) which 1s of type Ball as shown in (2). 


That(foa) = foa—(foae SE) (2) 


Where, SE is Ball structuring element of size 5 pixels and ° is the opening operation and That (foa) is the image 
after performing Top Hat operation. It is then followed by Bottom Hat operation. Bottom Hat operation 
involves subtracting the result of Top hat operation from the closing operation performed on it using a Ball 
structuring element (SE) as shown in (3). 


Bhat(X) = (That (foa) . SE) 7 That (foa)) (3) 


Where, Bnat(x) is the results of the Bottom Hat operation, SE is Ball structuring element of size 5 pixels and 
is the Closing operation. It is then contrast enhanced using (4). 


Ice = fod + That- Bhat (4) 


Where, Icg is the Contrast Enhanced image, That is the Top Hat Result and Bhat is the Bottom hat result. The 
contrast enhanced image is then median filtered (Imea) and is then subtracted from contrast enhanced image 
for complete background removal process using (5). 


Ign = linea - Ice (5) 


Where Ipn is the background removed image. The image is then contrast stretched to obtain the full dynamic 
range of the retinal image (Ics). It is then median filtered (Ima) again using a structuring element. H-maxima 
transform is then applied to the median filtered image (Ima) to suppress all the maxima in the intensity image 
Ima that is less than the threshold h, to obtain the image fp. It is then thresholded with the value of 0.05 
because features are not visible if the values are above 0.05. 


finx = Threshold(fp,0.05) (6) 
Where, finx is the thresholded image and fis the h-maxima Transform. Thus, the resulting image consists of 
Blood vessels, MA and HA. Blood vessels are detected using a modified method of Local entropy 
Thresholding and Morphological method. Firstly, local entropy Thresholding is used for the detection of the 
blood vessels, where the Optimal Threshold is calculated as shown in the following equation. 


Te= [maxrZg Hr(Th)] (7) 
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Where, Hr (Th) is the Second Order Entropy and Tx is the optimal threshold based on which the vessels are 
thresholded from the image. Detection of blood vessels (fpyi) 1s done using local entropy Thresholding 
method. Secondly, morphological detection of the blood vessels is done to enhance the detection of minor 
vessels. Morphological Blood Vessel detection involves the following steps. The green component of the 
image is extracted (f,). It 1s then inverted to obtain better enhancement. 


f, = 255-f, (8) 


Where, fj is the inverted image and f, is the green component image. After inversion of the image, canny edge 
detection is performed, which is a multistage algorithm to detect various edges. 


f. = canny(fi) (9) 


Where, f, indicates Canny detected image and f; is the inverted image. The image is dilated as well as eroded. 
The eroded image (f-) is subtracted from the dilated image (fq) to detect the border using the disk-shaped 
structuring element of radius 10 pixels. 


fo = fa-fe (10) 


Where, f, is the Border detected image, fa is the dilated image and f, is the eroded image. Contrast Limited 
Adaptive Histogram Equalization (CLAHE) is then applied on the green component image (fa). A 
morphological opening operation was performed on the image using structuring element (SE) of ball shape as 
given by (11). 


fo= fae SE (11) 


Where, f, is the result of opening operation, f.) is the concatenated CLAHE image and SE is the structuring 
element of disk shape. The opened image (f,.) then subtracted from the CLAHE (fc) performed image and 
then thresholded (fin). It is then median filtered for preserving the edges. The obtained image is then 
subtracted from the Border detected image (fp). The resultant image contains the detected Blood Vessels 
(fpbv2). Then, the results of both the blood vessel detection methods are combined to obtain enhanced blood 
vessel detection 
as given by (12). 


frinalbv = fovitfov2 (12) 


Where, ffinatby 1S the final detected blood vessel, fpyi is the blood vessel detected using local entropy 
thresholding and fpy2 is the blood vessel detected using morphology. The resultant blood vessels are again 
dilated using a disk SE. And finally it is subtracted from finx. 


frinal = ftnx-feinatbv (13) 


Where, frinatby 1S the detected blood vessel, fix is the image containing MA, HA, blood vessels and ffinai is the 
final MA, HA detected image. Thus, the image containing MA and HA are detected in the image. They are 
then separated based on the number of pixels into MA and HA. Hence, after the detection of MA and HA, 
features are extracted for classification. 


3.1. Feature Extraction 

Feature extraction helps in efficiently representing the interesting portions of an image as feature 
vector. Features like Area, Perimeter, Eccentricity (Ecc), Major Axis Length (MAL), Minor Axis Length 
(MIL), Energy (E), Contrast (C) and Homogeneity (H) are extracted. Area is the number of white pixels 
existing in the region of interest in the binary image. Perimeter is the number of white pixels present in the 
boundary of the region of interest in the binary image. 
Eccentricity is the measure of deviation of the conic section from being circular. It is given by, 


Ex =(1-™ (14) 


Where, m - semi-major axis and n-semi-minor axis and e is the eccentricity. 
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Contrast is the separation between the brightest and the darkest areas of the image. It is calculated as, 
C=y,¥, Di, j (15) 


Where, Dli,jl is the Normalized GLCM. 
Energy is the sum of the elements that are squared in the matrix. It is calculated as, 


E= Yidj( —f)*Dii, jl (16) 
Homogeneity measures the proximity of the element’s distribution in the matrix. It is given by, 


D|ij| 
1+|i-j| 





H = di dj (17) 


Along with these features, total number of MAs and HAs are calculated. 


3.2. Classification 

Neural Network is used for Classification after the extraction of the features. Quick propagation 
Neural Network is used. The numbers of hidden Layers used are two and activation function used is Sigmoid 
Function. The method for Linear and Multi class classification of early Diabetic Retinopathy using Alyuda 
Neuro Intelligence Software are discussed below. 

Firstly, the features extracted from the Training set data are stored in the Excel Sheets in Comma 
separated value (CSV) Format. One Excel Sheet contains the total number of Microaneurysms and 
Haemorrhages along with type of grading according to Table 1 as suggested by Dupas et al. [22]. Another 
Excel Sheet contains the six features extracted each for Microaneurysms and Haemorrhages. Secondly, after 
the creation of the Test data set using Excel sheet, Alyuda Neuro Intelligence Software is used. Thirdly, the 
test data file (Excel Sheet) containing the number of MAs and HAs along with the respective grading criteria 
specified in the last column, is loaded for Multi class classification of early DR 1.e., Normal, Mild, Moderate 
and Severe. The Test data file containing 12 features of MA and HA is loaded for Linear classification 
seperately i.e., Normal and Abnormal Early DR Images.Finally, the data is analyzed and pre-processed. It is 
then followed by selecting the design, where the type of neural network is chosen as Quick Propagation 
Neural Network, and the number of Hidden Layers as two for training. It is then tested to see the classifier 
results. And the test results displayed on the GUI is fed for linear and Multi-class classification respectively 
for the classification results to be displayed. 


Table 1. Rule Based Grading of Early Diabetic Retinopathy 


Early Diabetic Retinopathy Stadium Number of Microaneurysms Number of Haemorrhages 
Normal 0 0 
Mild 1<MAs<5 0 
Moderate 5<MAs<15 0<HAs<5 
Serious MAs=15 HAs>5 


3.3. Patient Database Management System 

Patient Database management system is created using MySQL for storing the details of patients. It 
has records related to Patient ID, Name, Age, Weight, Gender, and Date, Patient History like Diabetes, 
Glaucoma, Cataract and the Count of Microaneurysms and Haemorrhages obtained by the hybrid method. A 
Graphical User Interface is created using MATLAB R20l6a. It is built using push buttons, popup menus, 
panels, edit text boxes, static text boxes and axes. It consists of display options for the various stages of the 
algorithm for the detection of MA and HA along with the Execution Time, Classification using Neural 
Network and Patient Database Management System. 


4. RESULTS AND DISCUSSIONS 

The Retinal fundus image is obtained as shown in Figure 3(1), The Green component of the image is 
taken with the Optic disk masked image overlaid on it as shown in Figure 3(11), The Top and Bottom Hat 
Transforms are applied on the image for background removal. Firstly, Top Hat is performed and then 
Bottom-Hat is performed as shown in Figure 3(i11). Background removal along with contrast enhancement is 
done for eliminating uneven variations as shown in Figure 3(iv).To increase the dynamic range of the retinal 
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image, contrast stretching is done as shown in Figure 3(v) and it is then median filtered, thus preserving all 
the edges and reducing the salt and pepper noise as shown in Figure 3(v1). h-Maxima is then applied on the 
image to suppress the values less than a definite intensity as shown in Figure 3(vi1). It 1s then thresholded as 
shown in Figure 3(viii). The blood Vessel is detected using Local Entropy thresholding and morphological 
method as shown in Figure 3(ix). After the detection of blood vessels by the hybrid method, it is subtracted 
from Figure 3(viii) to obtain HAs and MAs. Number of pixels present in the connected regions distinguishes 
HA and MA. Finally, HAs and MAs are detected as shown in Figures 3(x) and 3(x1). 





ae J 


(viii) 


(xi) 





(x) 


Figure 3. (1) Retinal Fundus Image (11) Optic Disk mask overlaid on Green component image 
(111) After Top and Bottom Hat Transform (iv) Background removal (v) Contrast Stretching (vi) Median 
Filtered (vii) h-Maxima Transform (vii)Thresholded (ix) Blood Vessel detection based on Hybrid method 
consisting of local entropy thresholding and morphological method (x) Haemorrhages (x1) Microaneurysms 


4.1. Results of Classification using Neural Network 
4.1.1 Results of Linear Classification for Early Detection of Diabetic Retinopathy 

Various features like Area, Perimeter, Eccentricity, Major Axis, Minor axis, Contrast are extracted 
from Microaneurysms and Haemorrhages of the test data using Hybrid method for Early Detection of DR. 
The features extracted are entered in the Microsoft Excel sheet and the last column contains the classification 
of the training set as normal or abnormal DR for linear classification as shown in Figure 4. It is then loaded 
into Alyuda Neural Intelligence software where the neural network is analysed, pre-processed, designed and 
trained. The Quick propagation neural network is chosen, which uses sigmoidal activation function. The 
numbers of hidden layers used are two. It is then tested for 260 images, which are classified into Normal or 
Abnormal as shown in Figure 5. 


[yd ’ Book5.csv = Microsoft Excel non-commercial use Ss Ea 
M16 - fi 0 
A 8 c D E F G H ! Jj + L at ona N Oo i) a R s if UE 

1 Subject Ni Area Perimeter Eccentrici! Major axis Minor axicontrast Areal Perimeter: Eccentrici! Major Axi: Minor Axi Contrasti Classify 

2 subi 4 7.3 0.7 7.1 5.1 0.048 42 30.2 0.97 15.6 3.6 0.12 Abnormal 

3 sub2 0 0 0 0 0 0 0 0 0 0 0 0 Normal 

4 sub3 i 9 0.1 1.15 1.15 0.0002 0 0 0 0 0 0 Abnorma 1 

5 subs 4 15.5 0.2 5.6 5.6 0.02 63 26 0.65 9.5 8.5 0.07 Abnormal 

6 subs 16 11.619 Cc $5.02 4.16 0.03 39 21 0.7 8.77 5.9 0.0996 Abnormal 

7 sub6 22 6.9 0.6 7.2 54 3 38 25 0.92 “4 a 0.13 Abnormal 

8 sub7 0 0 0 0 0 t*] 0 0 o 0 0 0 Normal 

9 subs 15 7.5 0.3 5.15 4a 0.04 35 238 0.78 9.6 3.38 0.09 Abnormal 

10 subs 8 8.24 0.78 4.2 2.67 0,009 104 39 0.63 13 10.2 0.02 Abnormal 

11 sub1o 9 74 0.6 7A 6.2 0.05 30 2 0.96 16 5 0.13 Abnormal 

12 sublli 8 3.5 0.5 7.8 6 0,03 40 29 0.8 14 13 0.13 Abnormal 
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“16 | subiS 0 o ] 0 0 0 0 0 o 0 OIN ormeal 

17 subi6 22 19 0.9311 9.38 3.4 0.012 87 42.24 0.92 17.9 6.97 0.011 Abnormal 
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28 sub27 12 9.69 0.4 4 4 0,0088 37 20 0.68 8.22 5.9 0,008 Abnormal 

29 sub28 21 7.1 0.6 6.5 4 0.03 38 28 0.7 14 6 0.13 Abnormal 
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44> ht Books «= 4 > 





Figure 4. Features Extracted for Linear Classification 
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Figure 5. The Test Image classified as an Abnormal Image after entering the features 


4.1.2. Results of Multi Class Classification for Early Detection of Diabetic Retinopathy 

Number of Microaneurysms and Haemorrhages of the test data obtained from the Hybrid method are 
entered in the Microsoft Excel Sheet. They are used as the Test data as shown in Figure 6. It is then loaded 
into Alyuda Neural Intelligence software where, it is analysed, pre-processed, designed and trained. Once the 
system is trained, it 1s tested for 260 Images and graded into normal, mild, moderate and severe as shown in 
Figure 7. 


book3.csv = Microsoft Excel non-commercial use 
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Figure 6. Features extracted for Multi class classification/Rule based grading 
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Figure 7. Results of Multi class classification/Rule based Grading 


4.2. Results of Patient Database Management System 

Patient Details like Patient ID, Name, Age, Weight, Gender, Date, Patient History like Diabetes, 
Glaucoma, Cataract and Count of Microaneurysms and Haemorrhages obtained by the Hybrid method are 
stored in the patient database created using MYSQL 5.5 as shown in Figure 8. Database created in MySQL 
was named as ‘patient’ and table created with all the records was named as ‘Details’. 
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Figure 8. Early DR Patient Database created using MySQL 


4.3. Results of Graphical User Interface created for Early Detection of Diabetic Retinopathy 

A Graphical User Interface is created using MATLAB R201l6a. It consists of display options for 
three algorithms for the detection of MA and HA along with the Execution Time, Classification using Neural 
Network and Patient Database Management System. The features extracted for the hybrid algorithm are 
displayed in the GUI. The patient details can be entered in the GUI, which is stored in the database. Figure 9 
shows the GUI of the Early DR for a retinal fundus Image which is graded as moderate. It is observed that for 
the retinal image loaded, the number of MA’s are 7 and the number of HA’s are 3. 
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Figure 9. Graphical User Interface of the Retinal Fundus Image graded as moderate 


Table 2 shows the performance analysis of the hybrid approach and the results for different 
databases using three different algorithms are presented in Figure 10. 
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Figure 10. Performance Charts of the three different algorithms for different databases 


The results of the work were compared with the ground-truths available for each image of the 
DIARETDB and MESSIDOR databases and validation was carried out for the hospital images by a Vitreo- 
retinal specialist. The outcomes shows that the values of accuracy and precision have increased using the 
hybrid approach and for each fundus image the execution time is only 7-9s which is substantial. Table 3 
illustrates the comparison of the 3 methodologies for the Detection of red lesions.In the Morphological 
method, some of the unique features identified are the division of an image into four quadrants before 
performing the pre-processing along with the masking of the optic disk using the centroid method. In the 
Entropy based method, some of the unique features identified are local entropy thresholding performed for 
detecting and eliminating the blood vessels, along with the active contour method for the detection and 
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elimination of OD to detect the early DR lesions. In the Hybrid method, a unique combination of h-maxima 
transform, top-bottom hat transform, filling formed a part of pre-processing along with a modified 
combination of entropy thresholding and morphological method is used for the accurate detection of the 
blood vessels. Comparison of state-of-the art techniques for red lesion detection with the proposed method 1s 
presented in Table 4. 


Table 2. Analysis of the Hybrid Approach 
Total no. of fundus images DR Normal TP FP TN FN 
260 192 68 192 16 a2 0 


Table 3. Comparison of the Three Methodologies for the Detection of Red lesions 
Retinal Fundus Image Morphology Method(MM) Local Entropy Hybrid Method(HM) 
Thresholding(LET) 
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Blood Vessel Detection is Vessel Detection is Accuracy of Vessel Detection 
not accurate in comparison improved in comparison is far better when compared 
with the Ground Truth. with the Ground Truth. with the Ground Truth(than 
Minor Vessels are not Minor Vessels are detected. the 2 methods).All the blood 
detected. vessels are detected. 





Optic Disc detection is Optic Disc Detection 
Optic Disc detection is accurate (Centroid Method  (Centroid Method) is accurate. 
accurate (Centroid Method). and Active Contour 
Method). 





More number of False Reduced number of False Less number of False 


positives for positives for positives for Microaneurysm. 
Microaneurysms. Microaneurysm. 


Table 4. Comparison with Other State of the Art Techniques for Red Lesion Detection 


Authors Execution Time No. of images Accuracy 
Kande et al[109] 35s 89 93.23% 
Syna Sreng et al[84] 9.53s - 90% 
Hassan et al[54] lmin 60 89% 
Proposed method 7-9s 260 93.9% 


5. CONCLUSIONS 

A characterization system for red lesions to facilitate early DR detection was developed. 260 retinal 
images obtained from different databases and hospitals were considered for testing. The method achieved an 
accuracy and precision of 93.9% and 92.5%, respectively. The execution time was found to be 7-9 seconds. 
The automated system was developed using MATLAB R2016a and Alyuda Neuro Intelligence Tool was 
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used to classify the retinal fundus images into normal and abnormal and grading of the retinal images based 
on the presence of the MA and HA as normal, mild, moderate and severe. A patient database management 
system created using MySQL helped in storing all the patient details for future reference and a user friendly 
Graphical User Interface was developed for the Early Detection of Diabetic Retinopathy. 
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